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INTRODUCTION 


Itost  fluids  -iiian  cooled  within  a  certain  range  below  their  noraal  aeltlng 
point  will  reaiain  fluid  for  an  Indefinite  period  of  tlae  If  not  disturbed.  If 
8  sample  of  a  fluid  In  such  a  oetastable  condition  is  disturbed  by  the  intro<» 
ductlon  of  a  sead  crystal  It  will  solidify.    If  the  saapla  Is  snail,  such  as 
that  held  by  a  watch  glass,  the  solidification  or  crystallization  of  som 
fluids  Is  seeadagly  Instantaneous.    This  phenosenon  is  noted  for  benzene 
(CgHg)  as  well  as  for  many  other  substances. 

When  aueh  a  saii^>le  Is  undercooled  In  a  large  enough  quantity  and  distrlb- 
utad  as  in  a  long  tuba,  the  velocity  with  which  crystallization  proceeds  after 
**trlsgttriag"  is  observed  to  be  finite.    Fxirthermore  the  velocity  is  observed 
to  vary  directly  with  the  degree  of  undercooling. 

In  this  escperlment  an  attempt  is  made  to  show  how  the  velocity  of  crys* 
talllzatlon  of  benzene  varies  with  degree  of  undercooling.    The  results  will  be 
eotspared  to  a  similar  experlo^nt  performed  with  liquid  phosphorous  Powell, 
at  al,  (4)  and  with  theory.    An  empirical  equation  Is  then  proposed  idilch 
will  fit  the  data. 

K^UIFMSNT  AMD  PROCEDURS 

Tbe  saiqtle  holder  was  constructed  of  two  concentric  pyrex  glass  tidies  of 
l(km  and  6ibi  outside  diameters.    The  wall  of  the  outside  tube  was  approximate- 
ly 1mm  thick.    The  systcsi  wma  sealed  in  such  a  way  that  the  only  entry  was 
threuijb  the  side  arm  at  the  top.    Distilled  and  sodiias  dried  benzene  was  pi- 
petted through  this  side  arm  to  a  level  about  one  inch  above  the  bottom  end  of 
the  inside  tube.    This  arrangement  formed  at  the  top  of  the  aasiple,  a  one  inch 
iMll  lAiich  ceuld  be  used  to  initiate  freezing  without  contamination. 
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Ihft  cooling  bath  ctmtlAted  of  a  transparent  glass  container  filled  vith 
coniirelal  ethenol.    Ibe  container  had  to  be  constructed  in  the  oHumer  of  a 
Oewur  flask  in  order  to  proaMte  effici«at  cooling  and  to  prevent  frosting. 
The  outside  diameters  of  the  walls  of  the  container  vere  Slan  and  32nB.  Con- 
struction of  the  sss^le  holder  and  the  cooling  hath  is  shown  in  Plate  I. 

Xhe  hath  was  cooled  hy  passing  a  continuous  flow  ot  dry  air  through  a 
copper  coil  laowrsed  in  a  dry  ice-acetone  solution  and  th«i  throu^  a  copper 
U-tube  iaanrsed  in  the  ethanol.    The  teaperature  o£  the  bath  was  controlled 
by  varying  the  flow  of  air. 

Tesperature  observations  were  made  vith  a  double  junction,  eopper- 
e<»x8tantan  thenaocouple  and  a  Rt^icon  potentioneter.    One  Junction  of  tha 
thcxmcovple  was  enclosed  in  a  glass  sleeve  and  ioMrsed,  along  with  the 
•ai^le  holder,  in  the  ethanol  bath.    1%e  reference  junction  was  kept  in  a 
water-ice  bath.    A  stirring  aechanisB  was  devised  to  eliainate  the  tcaperature 
gradient.    The  pressure  on  the  saa^le  was  aaintained  at  760aa  by  connecting  a 
■snoewter  to  the  side  am  of  the  saatple  bolder. 

When  the  system  had  reached  an  equilibrium  teiq>erature,  freezing  was 
inltiatad.    This  was  done  by  introducing  a  small  piece  of  dry  ice  into  a  few 
drops  of  acetone  kept  in  the  bottom  of  the  well  of  the  sample  holder. 

Bat*  were  recorded  on  Dupont  1^,  931A,  Rapid  Reversal,  fine  grain,  hl^ 
spaad,  Mtion  picture  film.    The  camera,  a  Bell  and  Howell,  Filmo  70,  was 
placed  about  20  inches  from  the  location  of  the  saiple  holder  in  the  bath.  A 
millimeter  scale  was  taped  altmg  the  side  of  the  bath,  parallel  to  the  sample 
and  at  the  same  distance  from  the  camera.    The  camera  was  calibrated  with  a 
stroboscope  to  run  at  40  frasws  par  ••eond.    A  ten  second  sweep  clock  was 
included  in  the  picture  of  the  system  as  a  check  on  the  timing  for  all  runs. 
The  clock  revealed  that  the  camera  speed  tias  accurate  to  4;  1  frame  per  second 
throughout  the  experiment. 
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The  data  vere  taken  by  observing  the  film  with  a  w>vl*  filai  and  Klerofllm 
viewer  in  which  the  Individual  £r«M  vas  enlarged  to  apprwciaately  19  tiMM 
its  actual  size.    The  position  of  the  front  of  the  advancing  crystal  grovth 
was  recorcied  appreaciaately  4ca  on  either  side  of  the  location  of  the  thersM- 
couple  junction  in  order  to  ellainate  the  effect  of  any  temperature  gradient. 
The  tiae  of  this  transit  of  the  £ront  was  deterained  by  counting  the  nunber  of 
frsaes  Involved.    Velocity  was  determined  by  the  slope  of  a  plot  of  position 
versus  tia^.    In  seae  eases  a  slight  curvature  was  noticed,  which  was  probably 
caused  by  a  saall  teaverature  gradient.    In  any  case  the  velocity  was  deter- 
ained by  fitting  the  best  straight  line  to  the  points.    Hie  data  appear  in 
Sable  1. 

RESULTS 

The  velocity  of  cryatallizatlon  of  benzene  was  plotted  as  a  function  of 
teaqperature  as  shown  in  Plate  II.    Powell*  et  al.  (4)  obtained  a  siailar  plot 
for  liquid  phosphorous,  as  shewn  in  Plate  III.    The  velocity  does  not  appear 
to  approach  a  asxtaua  in  either  case,  for  the  tei^erature  ranges  used.  The 
results  for  liquid  phosphorous  were  closely  reproduced  by  the  authors  with  the 
equation  (V^)^/^  .  ^onst  (Tq-T)  where       Is  the  viscosity  and  T©  the  aelting 
point.    A  plot  of  (V-)^)^/2  versus  T  for  b«nB«M  slMrai  a  definite  curvature  as 
seen  in  Plate  IV.    This  is  cimtrary  to  the  results  for  liquid  i^osphorous. 

The  smomt  of  undercooling  achieved  by  Powell,  Gilaan  and  Hildebrand  for 
liquid  phosphorous  was  approxiaately  229  C  below  the  aslting  point  of  44. 1^  C. 
Below  that  teaperature  the  saaple  froze  spontaneously.    The  velocity  of  crys- 
tallization at  this  lowest  teaqperature  was  of  the  order  of  210  ca/sec.  Even 
at  this  high  velocity  the  authors  decided  that  the  sleeve  of  the  solid  next  to 


6 


Table  1.    TtteperaCures  In  ailllvolts  and  corresponding  temperature*  In  degrees 
centigrade*;  also  the  corresponding  velocities. 


Temperatiire 


Teaperature 

 i!S2  


Velocity 
(ca/sec) 


+0.102 

•f 0.025 

•0.040 

-0.096 

-0.150 

-0.204 

-0.166 

-0.251 

-0.226 

-0.188 

-0.131 

-0.114 

-0.070 

-0.008 

-I-0.041 

■I-0.086 

-0.216 

-0.245 

-0.172 

-0.087 

-0.052 

-0.026 

-0.003 

+0.015 

+0.055 

+0.075 


+2.6 

+1.5 

+0.6 

-1.1 

-2.5 

-3.9 

-5.3 

-4.3 

-6.6 

-5.9 

-5.0 

-3.4 

-3.0 

-1.8 

-0.2 

+1.0 

+2.2 

-5.6 

-6.4 

-4.5 

-2.3 

-1.4 

-0.7 

-0.1 

+0.3 

+1.4 

*1.9 


1.6 
1.6 
4.0 
5.7 
11.7 
15.6 
20.5 
17.8 
29.6 
28.0 
20.6 
14.0 
14.9 
8.4 
5.9 
2.6 
1.0 
23.4 
27.4 
15.6 
9.4 
5.9 
4.9 
5.3 
3.1 
2.3 
1.5 


*  Teaperatore  in  degrees  corresponding  to  that  In  millivolts  for  a  copper- 
constantan  thexvocouple  was  determined  from  data  by  Dike  (2) . 


ECPLANAXION  OF  PLATE  II 


Plot  of  velocity  of  crystallization  of  benzene 
▼ersiis  temperature. 

o    Data  taken  during  one  12  hour  period. 

X    Data  taken  3  days  later  using  tiMt  ■■■■ 
Mnple. 
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PLATE  II 


EXFLANAIIOM  OF  TUtXE  111 


Plot  of  vttlocity  of  crystallization  of  liquid  phof 
phorous  versus  teaperatura  frooi  work  by  Powell,  Gllaan 
and  Hlldebrand  (4). 
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T  in  ®C 


EXPLANATION  OF  PLATE  IV 

A  plot  of  (ytf)^^^  versus  taq>ecature  ^ere  th«  valu«« 
•f  V  ware  taken  from  a  bast  fit  of  the  experiaental  data  for 
benzene.    Values  of  >j  were  determined  from  an  extrapolation 
of  values  from  the  Handbook  of  Physics  and  Chemistry  (6) . 
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PLATE  IV 


IS 


the  bath  was  formed  at  nearly  the  bath  teoperature.    Their  decision  was  basftd 
on  the  ASsoBption  that  the  heat  o£  fusion  tor  the  solid  foraing  nearest  the 
bath  was  conducted  away  through  the  wall  o£  the  glass  tube  to  the  bath  rapidly 
MMi^  to  result  In  a  nearly  isothecaal  process.    The  front  of  crystallization 
near  the  center  of  the  tube  aust  advance  aore  slowly.    This  phenoacnon  is 
actually  observed  with  benzene.    Tho  above  assuaptimi  sesas  «wn  aore  probable 
in  the  present  case  since  the  aaxlana  velocity  attained  was  only  on  the  order 
of  30  ca/sec.    One  cannot  attribute  this  "sleeve"  pheoeaenon  to  a  teaperature 
gradient,  for  the  sleeve  was  (*served  for  approach  to  the  desired  teaperature 
either  frea  above  or  below  that  teaperature.    hoao  (I)  aentioned  in  his  work 
with  xylenols  that  the  velocity  of  crystallization  was  slightly  decreased 
v^m  Che  thickness  of  the  glass  tube  was  increased  froa  l.2Sm  to  a.OOaa, 
idiile  holding  all  other  paraasters  constant.    Considering  a  specific  exMple, 
the  velocity  was  decreased  by  about  7%  froa  3.6  ca/sec.    It  is  difficult  to 
iaagine  any  other  reason  for  this  than  that  the  conduction  of  the  heat  of 
fusion  frea  the  systea  is  aore  easily  acce^>lished  for  the  thinner  walls. 


DISCUSSION 


Tumbull  (5)  has  preeeatad  a  theory  of  phase  boundary  aigration  for  a 
one  coapeoeat  systea  where  the  effects  due  to  density  change  are  ignored. 
Consider  the  special  case  for  a  liquid  (a)  to  solid  (b)  transition  at  an 
uadorcooling  of  To-T.    A  aoleeule  in  order  to  cross  the  interface  froa  the 
liquid  to  the  solid  aust  acquire  an  activation  energy  Z^G.    The  xrequency  of 
transfer  of  a  aoleeule  in  this  direction  is 
(I)  (-^G/kT) 

lAere  j/©  is  the  frequency  of  the  asleoule  striking  the  "barrUr".    If  AG'  is 
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the  difference  In  Gibbs  potentials  per  iiolecule  between  liquid  and  solid,  the 

frequency  of  transfer  In  the  opposite  direction  Is 

(2)  l/ba  "Uo  exp[-(AG  +  nG*)/kT  ] 

The  net  frequency  of  transfer  of  laolecules  across  the  phase  boundary  Is 

<3)  }^  «Z4b  -X4a 

and  the  velocity  of  lalgratlon  of  this  boundary  Is  proportional  to  V.  Then 

(4)  V  -  C  exp  (- AG/kT)  [l  -  exp  (-AG7kT)j 
where  C  Is  the  constant  of  proportionality. 

From  the  basic  lavs  of  theraodynaalcs ,  If  we  neglect  density  changes  and 
consider  a  one  component  systcs  at  constant  pressure,  Glbbs  potential  can  b« 
expressed  as 

(5)  dGa  »  -S^dT  dG^  -  .S|,dT 

^ere  S^^  and  S],  are  the  entropies  referred  to  one  nwlecule  at  t«q>erature  T. 
Integrating  from  normal  equllibrluo  point  G(po,To)  to  G(po,T) 

(6)  G^-G®  -  -^MT  Gb-Gj  -  -  ^S^dT 
Subtracting  we  get                                        ^  * 

(7)  (G^-Gj,)  -  (gJ-gJ)  -  -  /(Sa-S^)  dT 
But  "^^ 

(«)  -  -  VT 

where  >  Is  the  latent  heat  of  fusion  per  molecule.    We  know  also  that  at  the 
aonsal  equlllbrluoi  point  G.  =  Gfa  so 

T 

(9)  G^-G^  -  AG'  -  AfdT/t 

and  finally 

(W)  AG'  -   Aln  T/To. 

The  resulting  expression  for  velocity  of  crystallization  In  an  undercooled 
liquid  Is  then 

(U)  V  -  C  exp  (ZM^       -  exp  JlAiSiTZTfti] 
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The  constants  C  and  AC  were  evaluated  froa  the  two  e<tuatlons  tAilch  result 
i^en  the  experlnental  values  of  V  at  tra&peratures  •A  and  -ffi  C  ware  8obstl<- 
tuted  in  equation  (11) .    It  was  found  that  the  expression  fits  the  experl^ 
■antal  data  reasonably  x«sll  for  C  «  exp  (-36.72)  and  AG  «  -l.SA  x  10'^^.  The 
expression  Is  plotted  on  Plate  V. 

A  negative  activation  energy,  however,  does  not  have  any  physical  steanlng 
whM  one  considers  its  definition.    Moreover,  the  value  for  C  is  unreasonable. 
It  is  therefore  difficult  to  give  the  theoretical  expression  any  physical 
interpretation.    At  best  equation  (11)  can  be  considered  an  eaplrlcal  equation 
which  approxlaiately  fits  the  data. 

The  fault  with  the  above  expression  may  lie  in  the  derivation  of  AG*. 
The  thenaodynaaics  of  reversible  processes  was  o^iloyed,  which  laplies  that  • 
naar  equillbrluB  condition  exists.    This  is  definitely  not  the  case  for  the 
solidification  of  an  undercooled  liquid,  which  say  account  for  the  unaccept- 
able values  of  C  and  AG  obtained  in  Turnbull's  equation. 

Prenkel  (3)  has  isq>lied  that  the  rate  of  a  change  in  phase  froa  liquid 
to  solid  antst  be  Inversely  proportional  to  the  viscosity.    This  was  born  out 
by  the  results  obtained  by  Powell,  Gllouin  and  Klldebrand  for  liquid  ^ospho- 
rous.    If  it  is  assuaed  that  the  velocity  of  crystallization  of  benzene 
*«P««ds  on  viscosity  in  this  same  way,  a  ■odlfication  of  Turnbull's  theory 
can  be  made.    This  can  be  done  by  replacing  the  quantity   Mn  T/Tq,  which  was 
derived  <«»  the  basis  of  reversible  theraodynaalcs ,  by  an  arbitrary  function  of 
the  temperature  i^lch  is  to  be  deteradned  froa  the  aeasured  velocities. 

If        is  assigned  the  value  of  the  activation  energy  associated  with 
diffusion,  Turnbull's  expression  can  be  written  «• 
(12)  V  »  ^  f(T) 


EXPLANATION  OF  PLATE  V 


Plot  of  V  versus  T  froa  th«  expression 

V  -  C  exp(-  ^G/Wt)      1-exp  "  ^^^^^^  | 

•uperiaposed  on  the  experimental  points  for  benzene  f^ere 

C  -  exp(-36.72) 
and  AG  -  -1.54x10'^^. 
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PLATE  V 


It 


^ere  £(I)  =  I  -  exp  j .    If  AG '  can  be  writtea  as 

(X3)  AG'  -  gCT^-T)  »  S»^^o-T>'^ 

»»»/ 

(lAere  the  sub  aust  begin  with  n  =  I  in  order  that  A  G '  vanish  at  I  =-  1^) , 
thm  a  Taylor's  ejtpansioa  of  f(T)  about  I„  can  be  written 

AG'  is  considered  as  the  positive  increase  of  the  barrier  AG  seen  by  a 
teolecule  in  the  solid,  and  aoreover,  it  oust  remain  positive  as  long  as  is 
greater  than  T.    As  T  approaches  vmxy  near  T<„  the  first  teem  in  the  series  •£ 
equation  (13)  becooes  the  doodnant  one.    Therefore  C]^  oust  be  a  positive  nuober. 
However,  if       is  positive  then  fron  equation  (14)  f(T)  and  therefore  V  wst 
become  negative  as  T  approaches  Tq.    This  is  not  a  possible  situation  so 
mist  be  zero  or  at  least  negligibly  saall.    Equation  (14)  then  becoMS 

(15)  f  (T)  -  ^  (To-T)2  .[g-j  .  ^]  (T,.T)3  *  . . . 

If  all  but  the  second  order  term  can  be  neglected  the  velocity  expression  is 
equivalent  to  that  found  for  liquid  phosphorous  above.    However,  it  has  been 
•hom  that  this  does  not  hold  for  benzene.    If  the  third  order  term  is  in* 
eluded  the  velocity  becooes  of  the  form 

(16)  V  -  £1  (B  AT^  +  /XT^) 

Evaluating  the  constants  C'  and  B  from  experimental  points  of  the  benzene  data 
V  becomes 

(17)  V  -  -^i^  (3.99  At2  +  AT^) 

7 

This  expression  fits  the  experimental  points  quite  well  as  shown  in  Plate  VI. 
Eqiiation  (17)  reaches  a  aaxiw  at  a  temperature  of  approximately  -90®  C. 
This  value  is  reasonable  since  there  is  no  indication  that  the  measured  veloc- 
ity is  approaching  a  «— yjuwim  at  -b^  C. 


EXPLANATION  OF  PLATE  VI 


Plot  of  V  versus  T  from  the  expression 
V  -  —  (B^T^  +  At3) 
•uperiaposed  on  the  experioental  points  £or  benzene  where 

C  -  0.0123 
and  B  -  3.99. 
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PLATE  VI 


n 


CONCLUSION 

It  would  sees  that  the  fon  of  Turnbull's  expression  is  correct  but  that 
the  derivation  of  AG'  by  aeans  e£  ordinary  theraodynamic  expressions  is  not. 
Ihe  assuatption  that  the  velocity  of  crystallization  is  inversely  proportional 
to  the  viscosity  of  the  liquid  appears  to  be  acceptable.  Equation  (17) » 
though  physically  reasonable,  aust  be  considered  eaqtirical  until  the  two  con- 
stants C  and  B  are  given  physical  aeaning. 

Special  note  should  be  taken  of  the  agrcenent  between  the  expression 
found  by  Gilaan,  Powell  and  Hildebrand  and  that  derived  in  this  paper  froa  the 
thMry  of  Tumbull.    The  only  difference  was  the  appearance  of  the  tera  AT^. 
A  correction    of  some  sort  might  be  expected  for  substances  \^ich  do  not  have 
the  spherical  synoetry  of  aonatotaic  aolecules.    It  would  be  of  Interest  to 
investigate  the  velocity  of  crystallization  of  an  undercooled  liquid  of  even 
aore  c<Mtpllcated  aolecular  structure. 


Thm  author  wlsh«s  to  express  his  sincere  appreciation  to  Dr.  Basil 
Cumutte  for  his  advice,  coxmsel  and  excellent  suggestions  throughout  the 
progress  of  this  research;  to  Dr.  I.  C.  Hisatsune  for  his  aumy  helpful 
suggestions;  and  to  other  asabers  of  the  staff  for  valuable  discussion. 
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It  WM  9k9wmi  that  undercooled  Imkummm  fr«MEM  At  a  fl&lta  rata.  It 
m»  aUe  obsarvaa  tbac  tha  rata  •£  fcaaelug  aapantia  la  a  diract        m  tiM 
4atraa  af  uadaraaollns. 

A  MHpia  a£  iMixatM  mu  cooled  to  varloua  tMparatwaa  balov  its  Mlting 
point  and  crystalllaatioa  mm  iultiatad.    Zha  progcaaa  of  tha  cryotal  growth 
m»  raaordod  vlth  a  aotiao  pletara  caaora  raaaiag  at  a  callbratad  apaad  of  40 
frama  par  aaaoad.    Xha  valoeity  of  tha  pragraaa*  aa  takaa  froa  tha  fila,  «m 
plettad  as  a  function  of  tha  tMparatura.    Xha  raatata  wura  comparad  with 
■iailar  work  dona  on  li<fttid  phoaphorous. 

A  traatamt  of  a  thaoratieal  axpraaaloa  baaad  aa  tha  thaory  of  raactloa 
ratas  raoultad  ia  aa  approocimta  fit  to  tha  aaparlMntal  data»  bat  gata  risa 
to  fiWittaiM  i4iich  eottld  not  ba  givan  a  aatiafaecory  interpratatloa .    A  mm- 
i^t  diffarant  traataont  of  too  thoory  racultad  ia  tha  aora  satisfactory 
•faatioai 


